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A NONLINEAR  FRACTURE  MECHANICS  APPROACH  TO  THE  GROWTH  OF  SMALL  CRACKS 
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SUMMARY 

Using  linear-elastic  and  nonlinear  fracture  mechanics,  many  investigators  have  tried 
to  explain  the  growth  of  small  cracks  in  plates  and  at  notches.  These  studies  have  con- 
centrated on  the  growth  of  small  cracks  ranging  in  length  from  10~2  to  1 mm.  On  the 
basis  of  linear-elastic  fracture  mechanics,  the  small  cracks  grew  much  faster  than  would 
be  predicted  from  large  crack  data.  Nonlinear  fracture  mechanics,  in  particular  the 
J-integral  concept,  and  an  empirical  length  parameter  have  been  used  to  correlate  small 
and  large  crack-growth  rate  data.  The  physical  interpretation  of  the  length  parameter, 
however,  is  unclear. 

Recently,  some  investigators  have  suggested  that  crack  closure  may  be  a major  factor 
in  causing  differences  between  growth  of  small  and  large  cracks.  The  purpose  of  the 
present  paper  is  to  use  an  analytical  model  of  crack  closure  to  study  the  crack  growth 
and  closure  behavior  of  small  cracks  in  plates  and  at  notches.  The  calculated  crack- 
opening stresses  for  small  and  large  cracks,  together  with  elastic  and  elastic-plastic 
fracture  mechanics  analyses,  are  used  to  correlate  crack-growth  rate  data.  At  equiva- 
lent elastic  stress-intensity  factor  levels,  calculations  predict  that  small  cracks  in 
plates  and  at  notches  should  grow  faster  than  large  cracks  because  the  applied  stress 
needed  to  open  a small  crack  is  less  than  that  needed  to  open  a large  crack.  These  pre- 
dictions agree  with  observed  trends  in  test  data.  The  calculations  from  the  model  also 
imply  that  many  of  the  stress-intensity  factor  thresholds  that  have  been  developed  in 
tests  with  large  cracks  and  with  load-reduction  schemes  do  not  apply  to  the  growth  of 
small  cracks. 

The  current  calculations  are  based  upon  continuum  mechanics  principles  and,  thus, 
some  crack  size  and  grain  structure  exist  where  the  underlying  fracture  mechanics  assump- 
tions become  invalid  because  of  material  inhomogeneity  (grains,  inclusions,  etc.). 
Admittedly,  much  more  effort  is  needed  to  develop  the  mechanics  of  a noncontinuum.  Never- 
theless, these  results  indicate  the  importance  of  crack  closure  in  predicting  the  growth 
of  small  cracks  from  large  crack  data. 
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NOMENCLATURE 

A^  coefficients  in  stress-intensity  factor  equation  (k  = 1,2) 

b location  of  concentrated  force,  m 

b^  dimensions  for  partially-loaded  crack  (k  = 1,2),  m 

c half  length  of  crack,  m 

c^  half  length  of  initial  crack  in  load- reduction  scheme,  m 

c half  length  of  crack-starter  notch,  m 

n 

d he  1 f length  of  crack  plus  tensile  plastic  zone,  m 

E Young's  modulus  of  elasticity,  MPa 

F boundary-correction  factor  on  stress-intensity  factor 

1/2 

K stress-intensity  factor,  MPa-m 

N number  of  cycles 

n number  of  bar  elements  in  crack-closure  model 

P applied  load,  N 

p applied  load  at  initiation  of  load-reduction  scheme,  N 

P maximum  apulied  load,  N 

crack-opening  load,  N 
R stress  ratio  (S  . /S  \ 

r radius  of  circular  hole,  m 


applied  stress,  MPa 


Sc  crack-closure  stress,  MPa 

applied  stress  at  initiation  of  load-reduction  scheme,  MPa 
a x,  maximum  applied  stress,  MPa 

lllaX 

Smin  minimum  applied  stress,  MPa 

SQ  crack-opening  stress,  MPa 

t specimen  thickness,  m 

V crack-surface  displacement,  m 

V displacement  of  crack-tip  element,  m 

w half  width  of  specimen,  m 

Wj  half  width  of  bar  element  at  point  j,  m 

x,y  Cartesian  coordinates 

x^  coordinate  location  of  element  i,  m 

a constraint  factor,  a = 1 for  plane  stress  and  a = 3 for  plane  strain 

1/2 

AK  stress-intensity  factor  range,  MPa-m  ' 

effective  stress-intensity  factor  range,  MPa-m1^2 

t.i.C“Zonc  corrected  cf  fcctivs  stress-"  in  tens  it.  j rector  rariye,  MPa—m 
AKfch  threshold  stress-intensity  factor  range,  MPa-m1^2 

ASeff  effective  stress  range,  MPa 

n material  constant,  n - 0 for  plane  stress  and  n «*  v for  plane  strain 

v Poisson's  ratio 

P length  of  tensile  plastic  zone,  m 

0^  stress  on  segment  of  crack  surface,  MPa 

oQ  flow  stress,  MPa 

ou  ultimate  tensile  strength,  MPa 

o yield  stress  (0.2  percent  offset),  MPa 

ys 

u length  of  compressive  plastic  zone,  ir. 


INTRODUCTION 

Most  experimental  and  analytical  studies  on  fatigue-crack  growth  have  been  conducted 
on  "large"  cracks  with  lengths  in  excess  of  2 mm.  However,  in  many  engineering  struc- 
tures, crack  growth  from  “small"  pre-existing  flaws  is  a major  portion  of  the  component's 
tatigue  life.  Numerous  investigators  (1-11)  have  observed  that  the  growth  characteristics 
of  small  fatigue  cracks  in  plates  and  at  notches  differ  from  those  of  large  cracks  in  the 
same  material.  These  studies  have  concentrated  on  the  qrowth  of  small  cracks  ranging  in 
length  from  10”2  to  1 mm.  On  the  basis  of  linear-elastic  fracture  mechanics  (LEFM) , the 
small  cracks  grew  much  faster  than  would  be  predicted  from  large  crack  data.  This  behav- 
ior is  illustrated  in  Figure  1,  where  the  crack-growth  rate  is  olotted  against  the 
linear-elastic  stress-intensity  factor  range,  AK.  The  solid  (sigmoidal)  curve  shows  a 
typical  result  for  a given  material  and  environment  under  constant-amplitude  loading. 

The  solid  curve  is  usually  obtained  from  tests  with  large  cracks.  At  low  growth  rates, 
the  threshold  stress-intensity  factor  range,  AKth,  is  usually  obtained  from  load- 
reduction  (K-decreasing)  tests.  Some  typical  expei imental  results  for  small  cracks  in 
plates  and  at  notches  are  shown  by  the  dashed  curves.  These  results  show  that  small 
ciacks  grow  faster  than  large  cracks  at  the  same  AK  level  ar.d  that  they  also  grow  at 
AK  levels  below  threshold. 

Using  nonlinear  fracture  mechanics,  some  investigators  have  tried  to  explain  the 
growth  of  small  cracks  in  plates  and  at  notches.  In  particular,  the  J-mtegral  concept 
and  an  empirical  length  parameter  [3,-i]  have  been  used  to  correlate  small-crack  and 
large-crack  growth  rate  data.  The  physical  interpretation  of  the  length  parameter,  how- 
ever, is  unclear.  More  recently,  several  investigators  [S-llj  have  suggested  that  crack 
closure  [12,13)  may  be  a major  factor  in  causing  the  differences  between  the  growth  of 
small  and  large  cracks. 
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Many  of  the  comparisons  between  growth  rates  for  small  and  large  cracks  have  been 
made  in  the  low  crack-growth  rate  region  associated  with  threshold  (AKt^)  for  large 
cracks.  Measurement  of  AKt^,  using  load-reduction  schemes  [14-17]  , establishes  a regime 
of  crack  sizes  and  applied  stress  levels  where  crack  growth  would  not  occur.  But  small 
cracks  and  cracks  from  pre-existing  flaws  have  been  observed  to  grow  [1-3,10]  at  stress- 
intensity  factor  levels  below  AKth,  as  shown  in  Figure  1.  Here  again,  several  investiga- 
tors [18-20]  have  experimentally  shown  that  crack  closure  is  causing  the  stress-intensity 
factor  threshold  under  load-reduction  schemes.  On  the  basis  of  these  results,  an  impor- 
tant question  concerning  the  growth  of  large  cracks  at  extremely  low  stress  levels  arises. 
Is  the  AK-rate  relationship  for  large  cracks  under  constant-amplitude  loading  given  by 
the  dash-dot  curve  in  Figure  1?  Admittedly,  a threshold  could  exist  for  large  cracks  due 
to  material  inhomogeneities  and  environmental  factors,  such  as  oxide  or  corrosion  products 
on  the  crack  surfaces  [21,22] . But  these  thresholds  must  be  obtained  in  tests  without  any 
load-interaction  effects. 
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To  adequately  resolve  the  differences  between  the  growth  of  small  and  large  cracks, 
crack -closure  effects  must  be  considered.  This  paper  is  concerned  with  the  analytical 
treatment  of  crack  closure.  A simple  strip-yield  model  of  crack  closure  was  developed  in 
Reference  23,  but  a more  useful  model  was  developed  in  References  24  and  25.  There  have 
also  been  several  other  attempts  to  develop  simple  analytical  models  of  crack  closure, 
but  a review  here  is  beyond  the  scope  of  the  present  paper  (see  Ref.  25) . 

The  purpose  of  the  present  paper  is  to  use  an  analytical  closure  model  [25]  to  study 
the  crack  growth  and  closure  behavior  of  small  cracks  in  plates  and  at  holes,  and  to  study 
the  closure  behavior  of  large  cracks  under  lead-reduction  schemes  used  in  threshold  test- 
ing. The  model  used  here  simulates  the  influence  of  plane-stress  and  plane-strain  condi- 
tions on  yielding  at  the  crack  tip.  The  model  is  based  on  the  Dugdale  model  [26]  , but 
modified  to  leave  plastically  deformed  material  along  the  crack  surfaces  as  the  crack 
grows.  Two  configurations  are  analyzed:  (1)  a center-crack  tension  specimen  and  (2)  two 

symmetric  cracks  emanating  from  a circular  hole.  A surface-crack  configuration  is  also 
considered,  but  It  is  analyzed  as  a center-crack  specimen  with  a modified  stress-intensity 
factor  solution.  The  crack-opening  stresses  for  these  configurations  are  calculated  as  a 
function  of  crack  length  for  constant-amplitude  loading  and  for  the  load-reduction  schemes 
used  for  threshold  testing.  An  improved  method  of  calculating  crack-opening  stresses, 
using  crack-surface  displacements,  is  also  presented.  The  calculated  crack-opening 
stresses  are  compared  with  experimental  measurements  made  under  a load-reduction  scheme. 

Because  specimens  with  small  cracks  are  usually  subjected  to  high  3tress  levels, 
Elber's  effective  stress-intensity  factor  range  [13]  is  modified  herein  for  plasticity. 
Crack-growth  rates  under  constant-amplitude  loading  are  correlated  with  the  plastic-zone 
corrected  effective  stress-intensity  f actor_range , AKef f,  over  a wide  range  of  stress 
ratios.  Using  the  crack-growth  rate  and  AKeff  correlations  for  large  cracks  in  several 
materials,  predicted  growth  rates  for  small  cracks  in  plates  and  at  circular  holes  are 
compared  with  experimental  rates. 


FATIGUE  CRACK-CLOSURE  ANALYSIS 

The  following  sections  include  a brief  description  of  the  analytical  crack-closure 
model  [25]  and  of  the  assumptions  made  in  the  application  of  the  model  to  the  growth  of 
small  cracks.  Those  sections  also  include  a description  of  the  plastic-zone  corrected 
effective  stress-intensity  factor  range  and  its  correlation  with  crack-growth  rates. 


Analytical  Crack-Closure  Model 

To  calculate  crack-closure  and  crack-opening  stress  during  fatigue  crack  growth,  the 
elastic-plastic  solution  for  stresses  and  displacements  in  a cracked  body  must  be  known. 
Because  there  are  nn  closed-form  solutions  to  elastic-plastic  cracked  bodies,  simple 
approximations  must  be  used  when  finite-element  solutions  [23]  are  inappropriate.  The 
Dugdale  modei  [26]  is  one  such  approximation.  The  crack-surface  displacements,  which  are 
used  to  calculate  contact  (or  closure)  stresses  under  cyclic  loading,  are  influenced  by 
plastic  yielding  at  the  crack  tip  and  residual  deformations  left  in  the  wake  of  the  grow- 
ing crack.  The  applied  stres.3  level  at  which  the  crack  surfaces  become  fully  open  (no 
surface  contact)  is  directly  related  to  contact  stresses. 

The  analytical  closure  model  developed  in  Reference  25  is  for  a central  crack  in  a 
finite-width  specimen  subjected  to  uniform  applied  stress  (Fig.  2(a)).  Ttiis  model  is 
extended  herein  to  cracks  emanating  from  a circular  hole  in  a finite-width  specimen  sub- 
jected to  uniform  applied  stress  (Fig.  2(b)).  The  model  is  based  on  the  Dugdale  model, 
but  modified  to  leave  plastically  deformed  material  in  the  wake  of  the  crack.  The  primary 
advantage  in  using  this  model  is  that  the  plastic-zone  size  and  crack-surface  displace- 
ments are  obtained  by  superposition  of  two  elastic  problems--a  crack  in  a plate  subjected 
to  a remote  uniform  stress,  S,  or  to  a uniform  stress,  o,  applied  over  a segment  of  the 
crack  surface.  The  stress-intensity  factor  and  crack-surface  displacement  equations  for 
these  loading  conditions  on  a plate  with  a central  crack  are  given  in  Reference  25.  The 
corresponding  equations  for  cracks  emanating  from  a circular  hole  are  given  in  Appendix  A. 

Figure  3 shows  a schematic  of  the  model  at  maximum  and  minimum  applied  stress.  The 
model  is  composed  of  three  regions:  (1)  a linear-elastic  region  containing  a circular 

hole  with  a fictitious  crack  of  half-length  c + ,,  (2)  a plastic  region  of  length  p. 


and  (3)  a residual  plastic  deformation  region  along  the  crack  surface.  The  physical 
crack  is  of  length  c - r,  where  r is  the  radius  of  the  circular  hole.  The  compressive 
plastic  zone  is  to.  Region  1 is  treated  as  an  elastic  continuum;  the  crack-surface  dis- 
placements under  various  loading  conditions  are  given  in  Reference  25  for  a central  crack 
and  in  Appendix  A for  cracks  emanating  from  a hole.  Regions  2 and  3 are  composed  of 
rigid-per fectly  plastic  (constant  stress)  bar  elements  with  a flow  stress,  o..  The  flow 
stress  is  taken  as  the  ultimate  tensile  strength  t/au)  herein  because  some  tests  were  con- 
ducted at  very  high  applied  stress  levels.  (If  a nonlinear  stress-strain  curve  had  been 
incorporated  into  the  model  [27]  , the  need  for  an  assumed  flow  stress  would  have  been 
eliminated.)  The  shaded  regions  in  Figures  3(a)  and  3(b)  indicate  material  that  is  in  a 
plastic  state.  At  any  applied  stress  level,  the  bar  elements  are  either  intact  (in  the 
plastic  zone)  or  broken  (residual  plastic  deformation) . The  broken  elements  carry  com- 
pressive loads  only,  and  then  only  if  they  are  in  contact.  The  elements  yield  in  compres- 
sion when  the  contact  stress  reaches  -0o.  Those  elements  that  are  not  in  contact  do  not 
affect  the  calculation  of  crack-surface  displacements.  To  account  for  the  effects  of 
state  of  stress  on  plastic-zone  size,  a constraint  factor  a was  used  to  elevate  the 
tensile  flow  stress  for  the  intact  elements  in  the  plastic  zone.  The  effective  flow 
stress  oa0  under  simulated  plane-stress  conditions  was  0o  and  under  simulated  plane- 
strain  conditions  was  3aQ.  The  constraint  factor  is  a lower  bound  for  plane  stress  and 
an  approximate  upper  bound  for  plane  strain  [25) . The  constraint  factor  was  assumed  to 
be  1,  1.73  (Irwin's  plane-strain  constraint  [28]),  or  3,  depending  upon  the  test  condi- 
tions. Plane-stress  conditions  were  assumed  for  high  applied  stress  levels  and  for  crack- 
growth  rates  higher  than  the  rate  at  transition  (flat  to  slant  crack  growth)  [29]  because 
the  plastic-zone  size  would  be  of  the  order  of  the  plate  thickness.  Plane-strain  condi- 
tions were  assumed  for  crack-growth  rates  lower  than  the  rate  at  transition  because  the 
plastic-zone  size  would  be  small  compared  to  plate  thickness.  Where  the  constraint  condi- 
tions were  unknown,  such  as  a small  surface  crack,  plane-stress  and  plane-strain  condi- 
tions were  both  tried. 

The  analytical  closure  model  was  used  to  calculate  crack-opening  stresses  as  a func- 
tion of  crack  length  and  lead  history.  The  applied  stress  level  at  which  the  crack  sur- 
faces are  fully  open  is  denoted  as  SQ,  the  crack-opening  stress.  These  stresses  have 
been  calculated  by  two  methods.  One  method  was  based  on  the  stress-intensity  factor  due 
to  contact  stresses  at  the  minimum  applied  stress;  the  other  was  based  on  crack-surface 
displacements.  Under  constant-amplitude  loading,  the  crack  closes  first  and  opens  last 
at  the  crack  tip,  so  the  stress-intensity  factor  method  was  used  to  calculate  SD  [24,25]. 
However,  under  some  variable-amplitude  loading  conditions  and  tne  load-reduction  schemes 
used  for  threshold  testing,  the  crack  closes  first  and  opens  last  away  from  the  tip. 

Under  these  conditions,  the  stress-intensity  factor  method  was  inadequate.  The  crack- 
surface  displacement  method  was  more  accurate  in  calculating  the  applied  stress  level  ac 
which  the  crack  surfaces  open  completely  under  general  cyclic  loading.  In  this  paper,  the 
displacement  method  was  used  to  calculate  SQ  and  is  presented  in  Appendix  B. 

The  closure  model  was  only  used  to  calculate  crack-opening  stresses.  In  turn,  the 
crack-opening  stress  was  used  to  calculate  the  plastic-zone  corrected  effective  stress- 
intensity  factor  range  and  then  the  crack-growth  rate. 


Plastic-Zone  Corrected  Effective  Stress-Intensity  Factor  Range 

Elber's  effective  stress-intensity  factor  range  [13]  is  based  on  linear-elastic 
stress-intensity  factors.  For  small  cracks  and  high  stress  levels,  the  plastic-zone  sizes 
are  no  longer  small  compared  to  crack  size,  and  linear-elastic  analyses  are  inadequate. 

To  correct  the  analyses  for  plasticity,  the  Dugdale  plastic-zone  length  (p)  is  added  to 
crack  length  (c) , like  Irwin's  plastic-zone  correction  [28].  Thus,  the  plastic-zone  cor- 
rected stress-intensity  factor  range  is 


(S 
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where  d = c + u and  p is  calculated  at  the  maximum  applied  stress,  Smax,  using  the 
equations  given  in  Appendix  C.  The  boundary-correction  factor  F accounts  for  the  influ- 
ence of  hole  radius  and  specimen  width,  and  is  given  by  the  product  of  the  functions  F^ 
and  f£.  These  functions  are  given  in  Appendix  A,  by  Eqs.  (6)  and  (22),  respectively. 
(This  parameter,  AK,  is  similar  to  Barenblatt's  cohesive  modulus  130]  under  cyclic 
loading . ) 

The  plastic-zone  corrected  effective  stress- in tens i ty  factor  range  was  obtained  from 
Eq.  (1)  by  replacing  Smin  by  SQ,  like  Elber's  modi f icat. on , and  is  given  by 
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Fatigue  Crack  Growth  Rates 

Instead  of  using  an  equation,  to  relate  the  crack-growth  rates  to  'K0 , , a table- 
lookup  procedure  was  used  in  this  paper.  The  primary  advantage  in  using  a table  is  that 
the  baseline  crack-growth  rate  da* c can  be  described  more  accurately  than  with  a 
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multi-parameter  equation,  especially  in  the  transitional  region  (flat  to  slant  crack  ^ 

growth).  To  illustrate  the  construction  and  use  of  the  table,  an  example  is  shown. 

Center-crack  tension  (CCT)  specimens  (w  = 80  mm)  were  used  to  obtain  crack-growth 
rate  data  on  2024-T3  Alclad  aluminum  alloy  sheet  (t  = 2 mm)  material  [5],  The  tensile 
strength  (ou)  and  the  assumed  flow  stress  (oc)  of  the  material  were  475  MPa.  The  maximum 
applied  stress  level  Smax  was  77.5  MPa  at  an  R ratio  of  0.01.  A plot  of  log  dc/dN 
against  log  AK  was  constructed.  The  AK  values  at  selected  crack-growth  rates  are 
listed  in  the  following  table: 


dc 

dN' 

mm/cycle 

AK, 

MPa -m1//2 

AK, 

MPa-m1//2 

s 

o 

S 

max 

AKeff. 

MPa-m1//2 

7.6E-06 

6.59 

6.70 

0.573 

2.90 

3.0E-05 

8.79 

8.95 

.596 

3.65 

4. 3E-05 

12.64 

12.86 

.597 

5.23 

6.8E-05 

14.29 

14.58 

.596 

5.95 

1.2E-04 

17.58 

18.06 

.600 

7.30 

3. 8E-04 

19.78 

20.20 

.599 

8.18 

1. 3E-03 

25.28 

25.88 

.594 

10.61 

7.6E-03 

36.27 

38.25 

.578 

16.30 

The  rates  were  selected  so  that  straight  lines  between  each  adjacent  data  point  would 
describe  the  AK-rate  data  using  a visual  fit. 

Next,  the  plastic-zone  corrected  stress-intensity  factor,  AK,  was  computed  from 
Eq.  (1)  with  the  plastic-zone  size  given  by  Eq.  (31).  The  constraint  factor  (a)  was 
assumed  to  be  unity  because  most  of  the  data  had  rates  higher  than  the  rate  at  transition 
(about  3 x 10"  mm/cycle) . For  small  AK  levels,  AK  is  given  approximately  by 


A closure  analysis  was  then  conducted  on  a CCT  specimen  subjected  to  a maximum 
applied  stress  ( Smax ) of  77.5  MPa  with  an  R ratio  of  0.01.  The  initial  crack  length 
was  1.5  mm.  Crack-opening  stresses  were  calculated  as  a function  of  crack  length  from 
the  model.  The  crack-opening  stress  for  each  value  of  AK  is  listed  in  the  table. 

Values  of  AKoff  were  computed  from  Eq.  (2).  The  table  of  dc/dN  against  AK?ff  forms 
the  crack-growth  rate  data  used  in  all  predictions  made  on  this  material  and  thickness. 
The  crack-growth  rates  at  intermediate  values  of  AKeff  were  obtained  by  assuming  a 

power-law  equation,  C|  AKe2f,  between  the  two  adjacent  data  points.  For  AKeff  values 
outside  of  the  table,  the  power-law  equation  between  the  two  closest  data  points  was 
extrapolated . 


APPLICATION  Q1'  THE  CRACK-CLOSURE  ANALYSIS 

The  analytical  closure  model  is  used  to  calculate  crack-opening  stresses  under 
constant-amplitude  loading  and  under  the  load-reduction  schemes  used  for  threshold  test- 
ing. Some  comparisons  are  made  between  calculated  and  experimental  c^ack-opening 
stresses  for  these  loading  conditions  on  a steel  and  an  aluminum  alloy. 

Crack-growth  rates  under  constant-amplitude  loading  for  large  cracks  (c  > 2 mm)  in  a 
steel,  an  aluminum  alloy,  and  a cast  (nickel-aluminum-bronze)  alloy  are  correlated  with 
the  plastic-zone  corrected  effective  stress- intensity  factor  range.  Usinq  these  correla- 
tions, the  growth  rates  for  small  cracks  in  plates  or  at  circular  holes  are  predicted  and 
compared  with  experimental  growth  rates. 


Constant-Amplitude  Loading 

To  show  how  the  crack-opening  stress  influences  the  local  crack-tip  displacements 
under  constant-amplitude  loading,  anu  consequently  the  crack-tip  damage  or  crack  growth, 
the  applied  stress  (S)  is  plotted  against  the  displacement  ( Vc  ) of  the  crack-tip  element 
in  Figure  4.  The  maximum  applied  stress  was  0.4’o  and  the  specimen  was  cycled  at  an 
R ratio  of  zero.  Recalling  that  the  elements  in  the  plastic  zone  are  rigid-perfectly 
plastic,  no  change  in  displacement  occurs  until  the  element  yields.  The  calculated 
crack-opening  stress  from  the  displacement  method  (Appendix  B)  is  shown  as  the  solid 
symbol  on  the  loading  branch.  These  results  demonstrate  the  significance  of  crack  clo- 
sure, in  that  ail  cyclic  loads  below  SQ  would  cause  no  plastic  deformation  at  the  crack 
tip  and,  presumably,  no  crack  growth.  These  results  also  show  that  the  crack  opens 
before  the  crack-tip  element  yields.  This  was  caused  by  the  finite  stress  concentration 
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in  the  mcdel  (elastic  stress  concentration  was  about  30) . If  the  model  had  had  an 
infinite  stress  concentration,  like  an  ideal  crack,  then  the  element  would  have  yielded 
when  the  crack  tip  opened.  Likewise,  compressive  yielding  occurred  only  after  some 
unloading.  The  crack-closure  stress  (Sc)  is  shown  as  the  solid  symbol  on  the  unloading 
branch.  Sc  is  the  stress  level  at  which  the  first  crack-surface  element  at  the  crack 
tip  closes.  In  this  example,  the  crack  was  not  allowed  to  grow  during  the  loading 
branch.  If  the  crack  had  been  allowed  to  grow,  then  closure  would  have  occurred  over  the 
newly  created  crack  surface  after  a smaller  amount  of  unloading. 


The  closure  model  was  exercised  under  simulated  plane-stress  and  plane-strain  condi- 
tions for  constant-amplitude  loading.  The  calculated  SQ  values  for  the  CCT  specimen 
remained  nearly  constant  over  a very  wide  range  of  crack  lengths  [25)  . This  constant 
value  is  called  the  "stabilized"  crack-opening  stress.  The  stabilized  crack-opening 
stresses,  normalized  by  Smax,  are  plotted  against  stress  ratio  for  various  stress  levels 
in  Figure  5.  The  calculations  were  made  on  CCT  specimens  made  of  an  aluminum  alloy,  but 
the  results  also  apply  for  steels  and  titanium  alloys  when  the  material  behavior  is 
elastic-perfectly  plastic.  At  any  R ratio,  the  S0/Smax  values  are  lower  for  higher 
values  of  a and  for  higher  values  of  stress  level  (Smax/0o).  The  influence  of  stress 
level  was  more  pronounced  under  plane-stress  conditions  (solid  curves)  than  under  plane- 
strain  conditions  (dashed  curves) . 


Threshold  Testing 

Because  many  of  the  comparisons  between  the  growth  of  small  and  large  cracks  have 
been  made  in  the  crack-growth  rate  region  associated  with  thresholds  for  large  cracks,  it 
is  important  to  know  whether  the  threshold  is  a material  behavior  or  is  caused  by  the 
load-reduction  scheme.  Some  load-reduction  schemes  used  to  obtain  AKt^  are  shown  in 
Figure  6.  The  ratio  of  current  load  (or  stress)  to  the  initial  load  (or  stress)  is 
plotted  against  crack  length.  The  crack  length  at  initiation  of  the  load-reduction  scheme 
was  20  mm  in  a CCT  specimen  (w  = 100  mm).  The  solid  curve  from  Bucci  [17]  (proposed  ASTM 
test  method)  is  based  on  a constant  rate  of  change  in  normalized  plastic-zone  size  with 
crack  extension  and  is  independent  of  the  R ratio.  The  normalized  K-gradient, 

(dK/dc)/K,  was  -0.08  mm~l,  as  recommended.  The  dashed  curves  show  the  load-reduction 
scheme  proposed  by  Ohta  et  al.  [19]  for  R = 0 and  0.8.  Their  scheme  is  based  on  a con- 
stant AK-gradient,  d(AK)/dc.  The  AK-gradient  ranged  from  -20  to  -100  MPa/mml/2.  The 
smallest  AK-gradient  was  used  in  Figure  6.  The  step  function  was  proposed  by  Robin  and 
Pluvinage  [16] . The  load  was  reduced  by  10  percent  after  0.2  mm  of  crack  extension.  The 
crack  extension  was  large  enough  for  the  crack  to  traverse  the  plastic  zone  created  by 
the  previous  load. 

As  previously  mentioned,  several  investigators  [18-20]  have  experimentally  shown 
that  the  stress-intensity  factor  threshold  under  load-reduction  schemes  can  be  explained 
by  the  crack-closure  behavior.  Some  typical  results  on  an  aluminum  alloy  are  shown  in 
Figure  7.  Minakawa  and  McEvily  [31]  conducted  a threshold  test  on  a compact  specimen  and 
measured  the  crack-opening  loads  as  the  AK  level  approached  AKth-  The  crack-opening 
loads  were  determined  from  a displacement  gage  at  the  edge  of  the  compact  specimen.  This 
location  was  remote  from  the  crack  tip.  For  high  AK  levels,  the  P0/Pmax  values  ranged 
from  0.15  to  0.35.  The  horizontal  line  is  the  calculated  Po/Pmax  ratio  from  the  closure 
model  under  constant-amplitude  loading  with  plane-strain  conditions  (o  = 3) . The  calcu- 
lated ratio  agreed  fairly  well  with  the  experimental  values.  As  AK  approached  AKth, 
the  Po/pmax  ratio  rapidly  rose  and  the  ratio  was  nearly  unity  at  threshold.  Thus,  the 
rise  in  crack-opening  load  explains  why  the  threshold  developed.  But  what  caused  the  rise 
in  crack-opening  loads?  A number  of  suggestions  have  been  advanced  to  explain  this  behav- 
ior. Among  these  are  the  mismatch  of  crack-surface  features  observed  by  Walker  and 
Beevers  [32]  in  a titanium  alloy;  the  corrosion  product  formation  on  the  crack  surfaces, 
as  observed  by  Paris  et  al.  [33]  and  Ritchie  et  al.  [34];  and  the  variation  in  the  mode 
of  crack  growth  with  stress-intensity  factor  level  as  reported  by  Ohtsuka  et  al.  [35], 
and  Minakawa  and  McEvily  [31].  The  mismatch  of  crack-surface  features  and  corrosion  pro- 
ducts on  the  crack  surfaces  could  cause  the  surfaces  to  come  into  contact  at  a higher 
load  than  the  load  for  a crack  without  the  mismatch  or  corrosion  products.  The  mode  of 
crack  growth  near  threshold  is  a combination  of  Mode  I and  Mode  II  (tensile  and  shear) 
instead  of  pure  Mode  I.  The  mixed-mode  crack  growth  causes  an  irregular  crack  surface 
profile  and,  consequently,  the  possibility  of  crack-surface  mismatch.  The  analytical 
treatment  of  crack  closure  due  to  crack-surface  mismatch  or  corrosion  products  on  the 
crack  surface  is  beyond  the  scope  of  the  present  paper. 

However,  now  consider  the  influence  of  the  load-reduction  scheme  on  crack  closure. 
Using  the  closure  model,  a simulated  threshold  test  was  conducted  on  a CCT  specimen  made 
of  steel  using  Bucci's  load-reduction  scheme  [17],  The  material  flow  stress  was  800  MPa 
and  the  constraint  factor  was  assumed  to  be  unity.  The  applied  load  (or  applied  stress) 
against  crack  length  is  shown  by  the  solid  curve  in  Figure  6 at  an  R ratio  of  zero. 

The  crack  length  at  initiation  of  the  load-reduction  scheme-  Cj  ! was  20  mm.  As  the 
applied  stress  was  reduced  with  crack  extension,  the  closure  model  predicted  a rise  in 
crack-opening  3tress.  To  find  out  why  the  crack-opening  stress  rose,  the  crack-tip  dis- 
placements during  a single  cycle  are  shown  in  Figure  8.  The  figure  shows  applied  stress 
normalized  by  the  maximum  applied  stress  plotted  against  displacement  of  the  crack-tip 
element,  similar  to  that  shown  in  Figure  4 for  constant-amplitude  loading.  Again,  the 
crack-opening  stress  is  shown  by  the  solid  symbol  cn  the  loading  branch.  Here,  the 
crack-opening  stress  is  higher  than  the  value  calculated  for  constant-amplitude  loading 
(shown  by  the  dashed  line).  In  contrast  to  the  closure  behavior  under  constant-amplitude 
loading,  in  the  threshold  test  the  crack  did  not  close  near  the  crack  tip  but  did  close 
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away  from  the  crack  tip,  as  illustrated  in  the  insert.  The  insert  shows  a schematic  of 
the  crack-surface  displacements  near  threshold.  The  shaded  region  around  the  crack 
depicts  the  residual  plastic  deformations  left  in  the  wake  of  the  growing  crack.  The 
largest  residual  deformations  occur  along  the  crack  at  about  cy  and,  consequently,  the 
crack-opening  stress  is  controlled  by  the  displacements  in  this  region.  These  displace- 
ments tend  to  prop  the  crack  open  at  the  tip.  Figure  8 shows  that  the  crack-tip  element 
does  not  change  in  displacement  until  after  the  crack  surfaces  are  fully  open.  During 
unloading,  S<1  shows  the  applied  stress  level  at  which  the  crack  surfaces  close  away 
from  the  tip  and  Sc  shows  the  applied  stress  level  at  which  the  crack  tip  closes. 

Ohta  et  al.  [19]  tested  CCT  specimens  made  of  HT80  steel  (t  = 10  mm)  over  a wide 
range  of  R ratios  (-1  to  0.8)  to  determine  crack-growth  rates  down  to  threshold  and  to 
measure  crack-opening  stresses.  They  used  the  Ohta-Sasaki  load-reduction  scheme  (dashed 
curves  shown  in  Fig.  6)  and  a small  crack-tip  displacement  gage  to  measure  SQ  values. 
Some  typical  results  on  crack-opening  stresses  for  R = 0 are  shown  in  Figure  9.  The 
So/smax  ratio  is  plotted  against  AK  level  during  the  load-reduction  scheme.  The 
crack-opening  ratio  was  constant  down  to  a AK  level  of  about  15  MPa-m1/2.  Near  this 
level,  the  ratio  rapidly  rose  as  the  applied  loads  were  reduced  and  a threshold  developed 
at  about  7 MPa-m^-/2.  The  solid  curves  show  calculations  from  the  closure  model  with 
a = 1 and  1.73.  The  plane-stress  calculations  are  in  good  agreement  with  the  experimental 
data  down  near  threshold.  Because  their  crack-closure  measurements  were  made  on  the  sur- 
face of  the  specimen,  their  results  should  be  close  to  those  for  plane-stress  conditions. 
Similar  comparisons  at  R = 0.8  and  -1  show  similar  agrtenh  nf  Li  rween  ' 1 inr -stress  calcu- 
lations and  measurements  (not  shown) . 

But  crack  closure  near  the  free  surfaces  of  a specimen  may  not  control  the  overall 
growth  rate  through  the  thickness  of  a specimen.  In  fact,  the  crack-opening  stresses  vary 
through  the  thiekj^ss  of  a specimen  [3f>  ] . Thus,  the  overall  growth  rate  at  low  IK 
values  may  be  controlled  by  a higher  constraint  factor  than  that  cf  plane  stress.  To 
determine  the  proper  constraint  factor,  a three-dimensional  analysis  would  be  required. 
However,  the  correlation  of  crack-growth  rate  data  at  various  stress  ratios  could  also  be 
used  to  experimentally  determine  an  approximate  constraint  factor  [37],  Using  a plot  of 
AK  and  rate  for  various  R ratios,  as  shown  in  Figure  10(a)  for  the  HT80  steel,  the 
spread  in  the  data  as  a function  of  R can  be  used  to  find  the  crack-opening  stresses 
needed  to  correlate  the  R-ratio  data  on  a AKejf  and  rate  plot.  Because  crack-opening 
stresses  are  a function  of  constraint  factor  (see  Fig.  5) , an  approximate  constraint 
factor  can  be  found  to  correlate  the  data.  For  the  10-mm-thick  steel  specimens  at  low 
AK  levels,  plane-strain  conditions  were  expected  to  exist  under  constant-amplitude  load- 
ing. Figure  10(b)  shows  a plot  of  AKeff  and  rate  for  the  same  HT80  steel  data  using  a 
closure  analysis  under  plane-strain  conditions  (a  = 3) . Curves  were  drawn  through  the 
data  for  each  R ratio.  Growth  rates  above  5 x 10-6  mm/cycle,  hasicallv  generated  under 
constant-amplitude  loading,  correlated  quite  well  under  plane-strain  conditions  even  for 
high  AK  values.  But  the  threshold  data  showed  systematic  variations  with  AKeff- 
Recalling  how  well  the  closure  analysis  under  plane-si tess  conditions  (a  = 1)  predicted 
threshold  behavior  under  a load-reduction  schei./a  (set  Fig.  9),  the  growth  rate  data  on  the 
HT80  steel  wore  replotted  in  Figure  11  using  a closure  analysis  with  a = 1 (Fig.  11(a)) 
and  a = 3 (same  as  Fig.  10(b)).  These  results  show  that  plane-stress  closure  is  con- 
trolling the  growth  rates  near  threshold  under  a load-reduction  scheme. 

Apparently*  the  plane -Stress  ’■egiC.ns  n»a*-  *h<>  free  surfaces  of  a specimen  f lay  an 
important  role  on  crack-closure  effects  under  load-reduction  senemes  and  variable- 
amplitude  loading.  McEvily  [38]  conducted  experiments  to  confirm  the  importance  of  crack 
closure  near  the  free  surfaces.  He  found,  in  a test  on  a 6061  aluminum  alloy  specimen 
(t  = 12.7  mm),  that  a spike  load  caused  significant  crack-growth  delay.  In  a second  test, 
he  carefully  machined  25  percent  of  the  thickness  from  each  surface  oi  the  specimen  after 
the  application  of  the  spike  load  and  found  very  little  crack-growth  delay.  Thus,  the 
ciack -Hiiosuie  tfttet  under  sp»ke  loading  is  predominantly  a surface  plrenbtfiTibni 

As  shown  by  Figures  7 through  11,  the  threshold  behavior  observed  in  load-reduction 
schemes  [15-17]  is  caused  by  a rise  in  the  crack-opening  stresses,  and  the  rise  is  caused 
by  residual  plastic  deformations  left  in  the  wake  of  the  growing  crack.  Thus,  load- 
reduction  schemes  create  aititically  high  values  ot  Stress- inteiis ity  lactot  thresholds , 
and  the  crack-growth  rate  data  generated  under  these  load-reduction  schemes  cannot  be 
used  in  an  LEFM  analysis  to  predict  the  growth  of  small  cracks.  To  obtain  crack-growth 
rate  data  at  low  AK  levels  using  large  cracks,  the  threshold  testing  procedures  must  be 
changed.  The  residual  deformations  on  the  crack  surfaces  near  the  beginning  of  a thresh- 
old test  must  be  eliminated  (see  insert  in  Fig.  8).  The  deformations  could  be  machined 
away  during  the  test  or  large  compressive  loads  could  be  applied  to  the  specimen  to  remove 
this  deformation.  The  specimen  could  also  have  been  side-grooved  to  eliminate  plane- 
stress  conditions  at  the  crack  front. 


Growth  ot  Small  Cracks 

In  the  following  sections,  the  crack-growth  rates  measured  in  specimens  with  large 
cracks  (c  > 2 mm)  are  correlated  with  the  piastic-zone  corrected  stress-intensity  factor 
range,  AK eff.  This  correlation,  a table  of  AKeff  against  dc/dN,  provides  the  crack- 
growth  rate  data  for  the  particular  material  and  thickness.  Using  the  closure  model  with 
an  assumed  constraint  factor,  the  crack-growth  rates  for  small  cracks  in  plates  and  from 
holes  are  predicted  and  compar  d with  experimental  rates.  All  initial  cracks  are  assumed 
to  have  no  previous  load  history  (no  residual  deformations) . 
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CSA  G40.ll  Steel.-  El  Haddad  [3]  tested  CCT  and  single-edge-crack  tension  (SECT) 
specimens  made  of  the  CSA  G40.ll  steel  under  cyclic  loading  at  an  R ratio  of  -1.  The 
material  yield  stress  was  370  MPa  and  the  ultimate  tensile  strength  was  510  MPa.  All 
specimens  were  2.54  mm  thick.  These  tests  were  conducted  to  study  the  growth  character- 
istics of  small  and  large  cracks. 

The  AK-rate  data  for  the  SECT  specimens  with  large  cracks  are  shown  in  Figure  12  as 
solid  symbols.  Because  the  closure  model  has  not  been  developed  for  SECT  specimens,  the 
data  were  assumed  to  be  equivalent  to  data  generated  on  CCT  specimens.  (El  Haddad's  d^ta 
on  CCT  specimens,  to  be  shown  later,  were  in  good  agreement  with  the  SECT  specimen  data.) 
The  solid  lines  through  the  SECT  specimen  data  in  Figure  12  show  the  baseline  (AK-rate) 
data.  To  determine  \Keff»  a closure  analysis  on  a CCT  specimen  (Smax  = 135  MPa;  a - 1) 
was  conducted.  The  corresponding  crack-growth  table  is: 


dc 

dN' 

mm/cycle 

AKeff' 

MPa-m1^2 

2. 24E-06 

6.53 

1.02E-05 

7.69 

4.70E-05 

9.78 

1.02E-04 

12.62 

1.93E-04 

17.29 

3.18E-04 

22.46 

6.60E-04 

34.20 

1.52E-03 

58.31 

This  table  was  used  in  all  predictions  made  on  the  CSA  G40.ll  steel. 

El  Haddad  [3]  tested  CCT  specimens  with  small  cracks  at  high  stress  levels  (207  and 
240  MPa).  The  crack-starter  notch  was  a very  small  hole  (r  = 0.19  mm),  but  no  details 
were  given  on  how  the  crack  was  initiated.  These  results  are  shown  as  the  open  symbols 
in  Figure  12  and  show  that  small  cracks  grow  much  faster  than  large  cracks  at  the  same 
AK  level. 

In  the  closure  analysis,  an  extremely  small  through  crack  (0.01  mm)  was  assumed  to 
exist  at  the  edge  of  the  starter  notch.  The  constraint  factor  was  assumed  to  be  unity 
because  the  specimens  were  subjected  to  high  stress  levels.  The  predicted  growth  rates 
are  shown  as  solid  curves  in  Figure  12.  The  predictions  are  in  fair  agreement  with  the 
experimental  results.  Both  experimental  and  predicted  results  show  a minimum  rate  at  a 
particular  AK  level.  In  the  analysis,  this  behavior  was  caused  by  the  transient  behav- 
ior of  the  crack-opening  stresses  shown  in  Figure  13.  Here  the  calculated  S0/Smax 
ratios  are  plotted  against  crack  length  for  the  two  tests.  The  calculated  crack-opening 
stress  was  initially  the  minimum  applied  stress,  but  the  opening  stress  rapidly  rose  and 
then  tended  to  level  off  as  the  crack  grew  under  cyclic  loading.  Truyens  (39)  has  experi- 
mentally shown  a similar  trend  in  SQ  with  crack  length  in  a ship  steel.  The  effective- 
stress  range  (ASeff),  indicated  in  the  figure,  was  Smax  - S0 • Crack-growth  rates  are 
directly  related  to  ASeff.  As  the  crack  grew,  the  decrease  in  ASeff  and  the  increase 
in  crack  length  (increase  in  stress - intens ity  factor)  combined  to  cause  the  unusual  behav- 
ior shown  in  Figure  12.  Truyens  (391  also  observed  a similar  trend  in  crack-growth  rates. 

El  Haddad  (3)  also  tested  specimens  with  cracks  emanating  from  larger  holes  than 
that  used  in  the  previous  example.  The  hole  radii  were  4.76  and  7.94  mm  in  70-mm-wide 
specimens.  The  maximum  applied  stress  was  135  MPa  at  an  R ratio  of  -1.  His  experimen- 
tal data  on  the  cracked-hole  specimens  are  shown  as  open  symbols  in  Figure  14  The  solid 
symbols  show  his  results  from  CCT  specimens  also  tested  at  an  applied  stress  level  of 
135  MPa.  These  results  show  tnat  small  cracks  rrom  holes  grow  much  faster  than  would  be 
predicted  from  using  I.EFM  analyses  with  large  crack  data. 

In  the  analysis,  the  baseline  (AK-rate)  data  were  obtained  from  large  cracks  in  CCT 
and  SECT  specimens.  The  solid  lines  through  the  CCT  specimen  data  in  Figure  14  show  the 
baseline  data.  Here  again,  no  details  on  the  crack-starter  notch  were  given  in  Refer- 
ence 3.  Therefore,  it  was  assumed  that  a very  small  through-crack  (0.01  mm)  was  at  the 
edge  of  the  holes.  The  solid  curves  in  Figure  14  show  predictions  f ron  the  closure 
analysis  with  i = 1.  The  predicted  rates  were  in  reasonable  agreement  with  the  experi- 
mental data.  For  cracks  emanating  from  holes,  the  transient  behavior  of  the  crack- 
opening  stresses  (similar  to  that  shown  in  Fig.  13)  and  the  plastic-zone  corrected 
stress-intensity  factor  range  combined  to  cause  the  trends  shown  m Figure  14. 

Alclad  2024-T3  Aluminum  Alloy.-  Broek  15]  tested  center-crack  and  cracked-h  le 
specimens  made  of  2024-T3  Alclad  aluminum  alloy  subjected  to  tensile  loadine  at  an  R 
ratio  of  0.01.  The  yield  stress  of  the  material  was  375  Mra  and  the  ultimate  tensile 
strength  was  475  MPa.  All  specimens  were  2 m thick.  The  ‘.K-rate  data  are  shown  in 
Figure  15  for  the  CCT  specimens  and  for  specimens  with  small  cracks  emanating  from  holes 
of  various  radii,  ranging  from  2.5  to  20  mm.  Again,  these  results  show  •'hat  small 
cracks  from  holes  grow  faster  than  would  be  predicted  using  .in  bid M analysis  with  large 
crack  data  on  CCT  specimens. 
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In  the  closure  analysis,  the  constraint  factor  was  assumed  to  be  unity  because  most  ' * J 

of  the  rates  were  above  the  rate  at  transition  (flat  to  slant  growth) . The  rate  at  tran-  / 

sition  was  estimated  to  be  about  3 x 10-3  mm/cycle  (first  knee  in  baseline  data  at  about 
9 MPa-ml/2) . The  development  of  the  crack-growth  rate  table  was  discussed  and  was 

included  in  the  section  on  "Fatigue  Crack  Growth  Rates."  The  predicted  results  for  small  \ 

cracks  at  holes  are  shown  as  solid  curves  in  Figure  15.  Solid  and  open  symbols  of  the 
same  type  denote  corresponding  predictions  and  experiments.  For  small  radii,  the  analysis 
underpredicted  the  growth  rates;  for  the  larger  radii,  the  analysis  overpredicted  the 
rates.  In  all  cases,  the  predicted  rates  were  within  about  a factor  of  2 of  the  experi- 
mental rates. 

Cast  Nickel-Aluminum-Bronze  Alloy.-  Taylor  and  Knott  [10,40]  tested  through-cracks 
and  small  surface  cracks  under  three-point  bending  at  an  R ratio  of  0.1.  The  material 
yield  stress  was  250  MPa.  The  ultimate  tensile  strength  was  not  reported  in  References  10 
or  40,  and  the  strength  was  estimated  to  be  550  MPa  [41].  Most  of  the  3urface-crack  tests 
were  conducted  with  a maximum  outer  fiber  stress  of  330  MPa  [40] . The  small  surface 
cracks  were  naturally-occurring  cracks  from  casting  defects  on  the  surface  of  the  speci- 
mens. Some  of  the  largest  defects  shown  on  photomicrographs  [10,40]  had  lengths  (2c) 
ranging  from  0.05  to  0.15  mm.  The  AK-rate  data  for  the  surface-crack  specimens  are  shown 
as  open  symbols  in  Figure  16.  Reference  10  did  not  give  the  particular  test  conditions 
for  each  surface-crack  test,  but  the  test  conditions  were  assumed  to  be  identical.  Thus, 
the  data  show  a lot  of  scatter  which  is  undoubtedly  due  to  the  cracks  growing  through  dif- 
ferent microstructure. 

The  solid  curve  showing  the  threshold  at  about  7.5  MPa-ml/2  was  obtained  from  bend 
specimens  with  large  through-cracks.  The  results  on  large  cracks  and  small  surface  cracks 
show  that  small  cracks  grow  at  AK  levels  well  below  the  threshold  established  from  large 
crack  tests.  (Again,  the  details  of  the  large  crack  tests  were  not  reported  in  Refer- 
ences 10  or  40  . ) 


Because  eL aok— y row th  rate  data  yoneidteu  under  load-reduction  schemes  (see  comments 
in  section  on  "Threshold  Testing")  cannot  be  applied  to  the  growth  of  small  cracks  [2], 
the  dash-dot  curve  was  assumed  to  be  the  baseline  data.  In  the  closure  analysis,  the 
constraint  factor  was  assumed  to  be  either  1 (plane  stress)  or  3 (plane  strain) . The 
crack-growth  rate  table  for  the  corresponding  constraint  factors  is: 


dc 

3n’ 

mm/cycle 

AKeff (a  = 1) , 
MFa-m1/2 

AKeff (a  = 3)  , 
1/2 

MPa-rn 

2. 54E-09 

1.62 

2.55 

2.77E-06 

5.88 

9.28 

6. 59E-06 

6.84 

10.78 

1.31E-05 

7.71 

12.18 

3. 19E-05 

9.23 

14.58 

7.59E-05 

11.07 

17.46 

The  closure  analysis  has  only  been  developed  for  through-the-thickness  crack  con- 
figurations. To  analyze  the  growth  of  a small  surface  crack  under  bending,  the  analysis 
for  the  CCT  specimen  was  used  with  modification.  In  the  analysis,  the  CCT  specimen  was 
subjected  to  an  applied  stress  level  equal  to  the  outer  fiber  bending  stress.  For  small 
surface  cracks,  this  assumption  is  adequate  [42).  The  crack-depth-to-crack-length  (a/c) 
ratios  for  the  experiments  were  reported  to  be  about  0.8  [10).  For  this  a/c  ratio  and 
for  small  surface  cracks  (a/t  nearly  zero),  the  average  st ress- i ntens i ty  factor  around 
the  crack  front  is  about 


AK  = A.S  . " c 0.7  (4) 

Thus,  the  stress-intensity  factor  for  the  CCT  specimen  was  multiplied  by  0.7  to  convert 
the  eqjation  to  a small  surface  crack  under  bending.  The  crack-opening  stresses  calcu- 
lated for  a center-crack  specimen  were  also  assumed  to  apply  for  the  surface-crack 
configuration. 

The  predicted  results  for  small  surface  cracks  are  shown  as  solid  curves  in  Figure  16 
for  the  respective  values  of  i.  The  plane-strain  predictions  were  closer  to  the  experi- 
mental results  rhan  the  plane-stress  predictions.  However,  under  plar.e-strain  conditions, 
the  initial  high  growth  rate  (solid  symbol)  was  considerably  lower  than  the  experimental 
rates.  This  discrepancy  may  be  due  to  the  crack  growing  in  some  "weak”  direction  in  the 
microstructure  m the  early  stages  [10]. 


CONCLUSIONS 

An  analytical  fatigue  crack-closure  model  was  used  to  study  the  crack  growth  and  clo 
sure  behavior  of  small  cracks  in  plates  and  at  circular  holes,  and  to  study  the  closure 


y _ / 0 behavior  of  large  cracks  under  load-reduction  schemes  used  in  threshold  testing.  Some 
comparisons  were  made  between  calculated  and  experimental  crack-opening  stresses  for 
these  loading  conditions  on  a steel  and  an  aluminum  alloy. 

Crack-growth  rates  unde*,  constant-amplitude  loading  for  large  cracks  (lengths  greater 
than  2 mm)  in  a steel,  an  aluminum  alloy,  and  a cast  nickel-aluminum-bronze  alloy  were 
correlated  with  the  plastic-zone  corrected  effective  stress-intensity  factor  range.  Using 
these  correlations,  the  growth  rates  for  small  cracks  in  plates  and  at  circular  holes  were 
predicted  and  compared  with  experimental  growth  rates. 

The  following  conclusions  were  obtained: 

1.  Under  load-reduction  schemes  used  in  threshold  testing,  the  rise  in  crack-opening 
stresses  near  threshold  is  caused  by  residual  plastic  deformations.  Thus,  load-reduction 
schemes  create  artificially  high  values  of  stress-intensity  factor  thresholds. 

2.  Near  threshold,  the  calculated  crack-opening  stresses  under  plane-stress  condi- 
tions ag-ec.  d well  with  experimental  measurements  on  steel  specimens. 

3.  Crack-qrowth  rate  data  generated  under  load-reduction  schemes  for  large  cracks 
cannot  be  used  with  linear-elastic  fracture  mechanics  analyses  to  accurately  predict  the 
growth  of  small  cracks. 

4.  Snail  cracks  (less  than  1 mm)  in  plates  and  at  holes  grow  faster  than  large 
cracks  (greater  than  2 mm)  at  equal  stress-intensity  factor  ranges  because  the  applied 
stress  needed  to  open  a small  crack  is  less  than  that  needed  to  open  a large  crack  and, 
consequently,  the  effective  stress  range  for  small  cracks  is  greater  than  that  for  large 
cracks. 

5.  Closure  analyses  based  on  plastic -zone  corrected  effective  stress-intensity  fac- 
tor ranges  and  crack-growth  rate  data  from  large  cracks  predicted  the  growth  of  small 
cracks  in  plates  and  at  holes  reasonably  well. 


APPENDIX  A— EQUATIONS  FOR  STRESS- INTENS ITY  FACTORS  AND  CRACK-SURFACE  DISPLACEMENTS 


The  analytical  closure  model  for  a crack  in  a finite  plate  [25]  and  for  cracks  ema- 
nating from  a circular  hole  in  a finite  plate  requires  the  stress-intensity  factor  and 
crack-surface  displacement  equations  for  the  two  elastic  configurations  shown  in  Fig- 
ure 17.  The  equations  for  stress-intensity  factors  and  crack-surface  displacements  for  a 
crack  subjected  to  various  loading  in  an  infinite  plate  (r  = 0)  were  obtained  from  the 
literature  [43).  These  equations  are  modified  herein  for  cracks  emanating  from  a circular 
hole  in  a finite  plate.  Some  of  the  approximate  equations  are  verified  with  boundary- 
collocation  [44,45]  and  conformal-mapping  [46,47]  analyses. 


Stress- Intensi ty  Factors  for  Cracks  Emanating  from  a Circular  Hole  in  an  Infinite  Plate 


Remote  Uniform  Stress.-  The  stress-intensity  factor  equation  for  the  configuration 
shown’ in  Figure  17(a),  with  w equal  to  infinity,  was  obtained  by  fitting  to  the  results 
from  References  45  and  47.  Ttie  equation  is 


K*  = KSFS  = S fid  F’f 
h ^ n h 


(5) 


where  K®  is  for  a crack  in  an  infinite  plate  without  a hole  and  F®  is  the  boundary- 
correction  factor  for  the  circular  hole.  The  equation  for  F^  is 


FS  . n _ r 

* h V 1 3 * n 


(6) 


where  n = 1 is  for  a single  crack  and  n = 2 is  for  two  symmetric  cracks.  The  fuuc- 


f.  0.707  * 0.765'-  - 0.282'“  * 0.74-3  * C.872'4 


(7) 


0. 358  - 


. . 4 2 5 '■  ‘ 


1.378 


5t  - 


(8) 


where 


= r c 


Fq . (5)  is  within  0.2  percent  o:  the  results  in  Refer- 


ences 45  and  41 


Concentrated  l rce.-  The  Stt cos-intensity  factors  for  two  cracks  emanating  from  a 

concentrated  torc-'s  in  an  infinite  plate, 


circular “hole  sub 


u . e 


d to  a sy.ttme! 


pair 
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Figure  18,  were  obtained  from  a boundary-collocation  analysis  by  Newman  (unpublished) . (j-  * t 

An  equation  was  fitted  to  the  results  and  is 


» 


*5  ■ « 


2?d 

^7rd(d2  - b2) 


(9) 


where  k£  is  for  a crack  in  an  infinite  plate  without  a hole  and  is  the  boundary- 

correction  factor  for  the  circular  hole.  The  equation  for  the  correction  factor  is 


F 


P 

h 


1 ♦ »i( 


(10) 


Ax  = -0.02X2  + 0.558X4 


(11) 


A2  = 0.221X2  + 0.046X4 


(12) 


where  > = b/d  and  X = r/d  for  X < y < 1 and  0 < X < 1.  Eq.  (9)  is  within  ±1  per- 
cent of  the  boundary-collocation  results. 


Partially-Loaded  Crack.-  The  stress-intensity  factors  for  the  configuration  and 
loading  shown  in  Figure  17 (b) , with  w - <»,  were  obtained  by  using  Eq.  (9)  as  a Green’s 
function  and  integrating  over  the  segment  b^  to  b2-  The  resulting  equation  is 


where 


Kh  = ir  ^ g<y.x) 


(13) 


G(Y,X) 


3A, 


2(1  - X)2J 


sin 


-1 


y 


b2/d 

b3/d 


(14) 


A^  and  A2  are  giver,  by  Eqs. 
same  form  as  Eqs.  (5)  and  (9), 


(11)  and  (12),  respectively. 
Eq.  (13)  is  rewritten  as 


To  get  Eq.  (13)  into  the 


(15) 


c 0 

where  K^,  is  for  a crack  in  an  infinite  plate  without  a hole  and  F^  is  the  circular 

hole  correction  factor.  The  function  F^  is  given  by 


F 


0 

h 


G(v,\) 


(16) 


where  G ( > , X ) is  given  by  Eq.  (14). 

Crack-Surface  Displacements  for  Cracks  Emanating  from  a Circular  Hole  in  an  Infinite  Plate 

The  exact  crack-surface  displacements  for  the  configurations  shown  in  Figure  17,  with 
w equal  to  infinity  and  hnJe  radius  equal  to  zero,  were  obtained  from  Westergaard  stress 
functions  given  in  Rc'erence  43.  These  equations  are  modified  herein  for  cracks  emanating 
from  a circular  hox<-.  "no  following  sections  given  approximate  displacement  equations  for 
remote  uniform  stress  and  for  the  partially-loaded  crack. 

Remote  Uniform  stress . - The  approximate  crack-surface  displacements  for  the  configur- 
ations shown  in  Figure  17(a),  with  w = •*•,  are 


V 


s 

h 


(17) 


s 

for  x £ d,  where  = v for  plane  strain  and  - = 0 for  plane  stress.  is  given 

by  Eq.  (fc).  In  the  region  near  the  crack  tip,  Eq.  (17)  is  very  accurate. 


KV> 


^ Partially-Loaded  Crack.-  The  approximate  crack-surface  displacements  for  the  config- 
uration and  loading  shown  in  Figure  17(b),  with  w = °°,  are 


vh  = K<x>  + VJ-*>:K 


(18) 


where 


£ 


V (x) 


2(1  - n 


1 - n.2)a 

ST (b  “ 


x)  cosh 


1 I d - bx 
d I b - x I 


>r72  2 . -1  /b  \ 

+ \d  - x sin  g-  1 


b = b2 
b = b. 


(19) 


for  [ x | < d,  and  f£  is  given  by  Eq.  (16).  Again,  the  equation  is  very  accurate  in  the 
region  near  the  crack  tip. 

Finite-Width  Corrections 

The  equations  given  in  the  preceding  sections  for  stress-intensity  factors  and  crack- 
surface  displacements  (Eqs.  (5)  to  (19))  are  for  cracks  emanating  from  a circular  hole  in 
an  infinite  plate.  But  these  quantities  are  influenced  by  the  finite  width  of  the  plate. 
Therefore,  some  approximate  finite-width  corrections  are  developed  herein.  The  stress- 
intensity  factor  for  a crack  in  a finite-width  plate  is 


K = K F 
«®  w 


(20) 


where  is  the  stress-intensity  factor  for  cracks  emanating  from  a circular  hole  in  an 

infinite  plate  and  Fw  is  the  finite-width  correction  for  the  particular  loading  condi- 
tion. Noting  that  the  elastic  crack-surface  displacements  in  the  region  of  the  crack  tip 
are  directly  related  to  the  stress-ir.tensity  factor,  it  is  proposed  that  the  same  correc- 
tion factor  be  used  for  displacements: 


V 


V F 

00  w 


(21) 


where  VM  is  the  displacement  for  cracks  emanating  from  a circular  hole  in  an  infinite 
plate,  again  subjected  to  the  particular  loading  condition.  Eq.  (21)  gives  very  accurate 
crack-surface  displacements  in  the  region  near  the  crack  tip. 

Remote  Uniform  Stress.-  The  approximate  boundary-correction  factor  for  two  symmet- 
ric cracks  emanating  from  a circular  hole  in  a finite-width  plate  subjected  to  uniform 
stress  [48],  <-s  shown  in  Figure  17(a),  is 


/ wr  \ 

scc 


( ’'d  \ 
' 2w  / 


(22) 


for  r/w  < 0.5  and  d/w  < 0.7.  Eq.  (22)  is  within  *2  percent  of  boundary-collocation 
results  [45].  The  equation  accounts  for  the  influence  of  width  on  stress  concentrations 
at  the  edge  of  the  hole  [49]  and  the  influence  of  width  on  stress-intensity  factors.  The 
crack-surface  displacements  given  by  Eqs.  (17),  (21),  and  (22)  for  r = 0 were  compared 
with  results  from  Reference  44,  which  used  the  boundary-collocation  method.  The  displace 
ments  were  within  2 percent  of  each  other  for  any  value  of  x for  d/w  < 0.7.  The  dis- 
placements near  the  crack  tip  were  very  accurate,  as  expected. 


Partially-Loaded  Crack.-  The  approximate  boundary-correction  factor  for  two  symmetric 
cracks  emanating  from  a circular  hole  in  a finite-width  plate  subjected  to  partial  loading 
on  the  crack  surface  (Fig.  17(b))  was  obtained  from  the  infinite  periodic  array  of  cracks 
(r  = 0)  solution  [43].  The  correction  factor  is  modified  herein,  as  was  done  in  Refer- 
ence 25,  and  is  given  by 


(23) 


3^  = sin  ( "b^/2w^/sin  (-d/2w) 


(24) 


*»  £_' 


APPENDIX  B- -CALCULATION  OF  CRACK-OPENING  STRESS 

The  following  section  describes  the  method  used  to  calculate  the  crack-opening 
stresses  in  the  analytical  closure  model.  To  find  the  applied  stress  level  needed  to 
open  the  crack  surface  at  any  point,  the  displacement  at  that  point  due  to  an  applied 
stress  increment  (SQ  - Sm^n)  is  set  equal  to  the  displacement  at  that  point  due  to  tie 
contact  stresses  at  Smyn.  This  calculation  is  made  at  the  centroid  of  all  element.'  in 
contact  along  the  crack  surface.  The  maximum  value  of  (SD  - Smin)  gives  the  applied 
stress  level  at  which  the  last  element  (in  contact)  separates  (or  opens) . This  approach 
is  possible  because  the  stresses  on  the  elements  are  elastic  until  the  crack-tip  element 
yields.  Crack-tip  element  yielding  occurs  only  after  the  crack  surfaces  open  (see 
Figs.  4 and  8) . 

The  displacement  at  point  i,  along  the  crack  surface,  due  to  an  applied  stress 
increment  is 

V?  = (S  - S . )f(x.  (25) 

i vo  min/  V i/ 

where 

f(*i) = 2^i n2)  f2  - *i  FhFw  (26) 

F®  and  F^  are  given  by  Eqs.  (6)  and  (22),  respectively,  with  d replaced  by  c.  The 
displacement  at  point  i due  to  contact  stresses  is 


where 


;i  = I v(vxj) 

j=n 


g(vxj)  = G(xi'xj) + G(-vxj) 


G(xi'xj) = -Ve---  (b2  - xi) cosh_1  ( 


c - b2xi 
c|b2  - xii 


- - x^)  cosh 


31  - Xi 


1 1 i \ C / \ C ' J ii  w 


bj  = x-j  - w;  and  bi  = Xj  + Wj  (Wj  is  the  half  width  of  bar  element  at  point  j).  (In 
the  closure  model,  10  elements  were  used  in  the  plastic  zone  and  elements  11  to  n were 
residual  plastic  deformation  elements  along  the  crack  surface  (see  Ref.  25  for  details).) 
F^  and  f£  are  given  by  Eqs.  (16)  and  (23),  respectively,  again  with  d replaced  by  c. 
Equating  Eqs.  (25)  and  (27),  and  solving  for  ( SD > , gives 


(Sc).  = Smin  + b Cjg(Xi'Xj)/f‘Xi 


for  i = II  to  n.  The  maximum  value  of  ( SQ  ).  gives  the  crack-opening  stress,  Sc 


APPENDIX  C — CALCULATION  OF  PLASTIC- ZONE  SIZE 
Center-Crack  Tension  Specimen 

The  piastic-zone  size  (t)  for  a crack  in  a finite-width  specimen  (Fig.  2(a))  was 
obtained  non  Rice  [50]  and  rederived  by  Newman  [25]  as 


p = 


sin 


1 


(31) 


For  o=l,  the  equation  reduces  to  the  expression  derived  by  Rice  from  an  infinite  peri- 
odic array  of  cracks  solution.  Newman  derived  the  equation  from  approximate  stress- 
intensity  factor  equations  for  finite-width  specimens. 


Cracked-Hole  Specimen 

The  plastic-zone  size  for  cracks  emanating  from  a circular  hole  in  a finite-width 
specimen  (Fig.  2(b))  was  determined  by  requiring  that  the  finiteness  condition  of 
Dugdale  [26]  be  satisfied.  This  condition  states  that  the  stress-intensity  factor  at  the 
tip  of  the  plastic  zone  is  zero  and  is  given  by  r 


K?  + K°  = 0 
n n 


(32) 


where 


and 


K?  = S 
h max 


/^d  F?FS 
h w 


(33) 
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F®  and  F^  are  givan  by  Eqs.  (6)  and  (22),  respectively.  Fft 
from  Eqs.  (16)  and  (23) , respectively,  by  setting  by  = c and 
Eqs.  (33)  and  (34)  into  Eq.  (32)  and  simplifying  gives 


and  were  obtained 

b2  = d.  Substituting 


.«?  F 
max 


F®FS  - no  Cl  • --  sin-1  F°F°  = 0 

n w o l it  V d > J h w 


(35) 


For  a given  applied  stress  Smax<  constraint  factor  a,  flow  stress  c;0,  specimen  width 
w,  hole  radius  r,  and  crack  length  c,  Eq.  (35)  was  solved  for  d by  using  an  iterative 
technique.  Noting  that  p is  equal  to  d - c gives  the  plastic-zone  size. 

The  plastic-zone  sizes  for  cracks  emanating  from  a circular  hole  in  an  infinite  plate 
subjected  to  remote  uniform  stress  are  shown  in  Figure  19.  The  plastic-zone  size  normal- 
ized by  c is  plotted  against  S/(ac0).  The  dashed  curve  (r  = 0)  shows  the  Dugdale 
mcdel  [26].  The  solid  curves  show  the  normalized  plastic-zone  sizes  for  various  c/r 
ratios.  For  c/r  equal  to  unity,  yielding  does  not  occur  for  S/(ao0)  ratios  less  than 
one-third  or  1/KT,  where  KT  is  the  stress  concentration  for  the  circular  hole. 
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Figure  11.-  Crack-growth  rates  as  a function  of  stress-intensity  factor  range 
for  various  R ratios  using  closure  analyses  for  simulated  plane-stress 
and  plane-strain  conditions. 


Figure  12.-  Comparison  of  experimental  and  predicted  crack -growth  rates  for 
small  cracks  in  center-crack  tension  specimens  subjected  to  high  stress 
love  1 s . 


Figure  15.-  Comparison  of  experimental  and  predicted  crack-growth  rates  for 
small  cracks  emanating  from  a circular  hole  in  aluminum  specimens. 


Figure  16.-  Comparison  of  experimental  and  predict'd  crack-urowfh  rates  for 
snail  surface  cracxs  in  specimens  subjected  to  a high  stress  1 *v«.  1 . 


